Intermittency in Dynamics of Two-Dimensional Vortex-like Defects 
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We examine high-order dynamical correlations of defects (vortices, disclinations etc) in thin films 
starting from the Langevin equation for the defect motion. We demonstrate that dynamical correla- 
tion functions F2n of vorticity and disclinicity behave as F2n ~ y^ /r'^" where r is the characteristic 
scale and y is the renormalized fugacity. As a consequence, below the Berezinskii-Kosterlitz-Thouless 
transition temperature F2n are characterized by anomalous scaling exponents. The behavior strongly 
differs from the normal law F2n ~ F2 occurring for simultaneous correlation functions, the non- 
simultaneous correlation functions appear to be much larger. The phenomenon resembles intermit- 
tency in turbulence. 
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INTRODUCTION 

It is well known that defects like quantum vortices, 
spin vortices, dislocations and disclinations play an es- 
sential role in physics of low-temperature phases of thin 
films. Berezinskii ^ and then Kosterlitz and Thouless 
recognized that there is a class of phase transitions in 
2d systems related to the defects. The main idea of their 
approach is that in 2d the defects can be treated as point 
objects interacting like charged particles. It is usually 
called Coulomb gas analogy. The low-temperature phase 
corresponds to a fluid constituted of bound uncharged 
defect- antidefect pairs, which is an insulator, whereas 
the high-temperature phase contains free charged par- 
ticles and can be treated as plasma. Correspondingly, in 
the low-temperature phase the correlation length is infi- 
nite whereas in the high-temperature phase it is finite. 
A huge number of works is devoted to different aspects 
of the problem, see, e.g., the surveys [!-§. The scheme 
proposed by Kosterlitz and Thouless can be applied to 
superfluid and hexatic films and planar 2d magnetics. 
It admits a generalization for crystalline films, see Refs. 

. There are also applications to superconductive ma- 
terials, especially to high-Tc superconductors, see, e.g., 
Ref. 0. 

The dynamics of the films in the presence of the de- 
fects was considered in the papers |l^Jl^ . In the works a 
complete set of equations is formulated describing both 
motion of the defects and hydrodynamic degrees of free- 
dom. Then, to obtain macroscopic dynamic equations, 
an averaging over an intermediate scale was performed. 
At the procedure the "current density" related to the de- 
fects was substituted by an expression proportional to 
the average "electric field" and to gradients of the tem- 
perature and of the chemical potential. The resulting 
equations perfectly correspond to the problems solved in 



the works [|Tl|,|2|. Unfortunately, at the procedure an 
information concerning high-order correlations of the de- 
fect motion is lost. That is the motivation for the present 
work where these high-order correlations are examined. 

We start from the same "microscopical" equations of 
the defect dynamics as was accepted in Ref. Q. Fol- 
lowing the works we focus mainly on the case when the 
motion of the defects is determined by the Langevin equa- 
tion describing an interplay between the Coulomb inter- 
action and the thermal noise. We believe that the ap- 
proach is correct for hexatic films (membranes, Langmuir 
films, freely suspended films). The situation is a bit more 
complicated for the vortices in superfluid fllms because 
of the Magnus force. Nevertheless, the equation for the 
vortices is close to the Langevin equation, see Ref. [ pT| . 
Similar equations can be formulated for the dislocations 
in crystalline films, see Ref. for the vortices in super- 
conductors in some interval of scales, see, e.g., Ref. 
and for the spin vortices in planar 2d magnetics. We 
will not consider the last cases here, though our scheme 
is, generally, applicable to the systems. Treating non- 
simultaneous correlation functions related to the defects 
one should take into account creation and annihilation 
processes also. For the purpose we use the Doi technique 
]l3t who demonstrated that dynamics of classical parti- 
cles involved into chemical reactions can be examined in 
terms of the creation and annihilation operators, like in 
the quantum field theory. 

We consider correlation functions of the "charge 
density" p (vorticity, disclinicity etc) provided that the 
so-called renormalized fugacity y is small. The inequality 
?/ <C 1 is satisfied for large scales in the low-temperature 
phase and probably in some region of scales above T^. 
In statics, the normal estimate ~ F2 is valid at the 
condition. Surprisingly, the non-simultaneous high-order 
correlation functions F2n appear to be much larger than 
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their normal estimate i^^. In the low-temperature phase 
the phenomenon reveals an anomalous scaling on large 
scales. The reason for such unusual behavior is that the 
main contribution to high-order non-simultaneous cor- 
relation functions is associated with rare single defect- 
antidefect pairs. The situation resembles the intermit- 
tency phenomenon in turbulence, see, e.g., Ref. It 
can also be compared with non-trivial tails of probability 
distribution functions in the physics of disordered mate- 
rials, see, e.g., Refs. |l5| , p^ . Some preliminary results 



were published in the paper |17|. 

Let us give a qualitative explanation of the phe- 
nomenon. To obtain a non-zero contribution to the corre- 
lation function F2„(ii, . . . , ri, . . . , r„) one must con- 
sider trajectories of the particles passing through the 
points ri, . . . ,r„ at the time moments ti, . . . The 
situation is illustrated in Fig. [|. The "single-pair" con- 
tribution has to be compared with a "normal" contribu- 
tion associated with a number of defect-antidefect pairs. 
Though the normal contribution contains an additional 
large entropy factor it has also an additional small fac- 
tor related to a small probability to observe a defect- 
antidefect pair with a separation larger than the core 
radius. As a result of the competition, the normal con- 
tribution appears to be smaller. To avoid a misunder- 
standing, let us stress that the arguments do not work 
for the simultaneous correlation functions. The reason 
is that trajectories of two defects cannot pass through 
n > 2 points simultaneously, see Fig. ^j. 

Our paper is organized as follows. In Section |l| we 
remind some basic facts concerning static properties of 
the 2d defects and their dynamics and then we shortly 
review the Doi technique [ p^ suitable for our problem. 
In Section ^ we develop a diagrammatic representation 
for dynamical objects and examine the two-particle con- 
ditional probability which is extensively exploited in the 
subsequent consideration. In Section |3| we demonstrate 
how renornialization of different parameters can be ob- 
tained in the framework of our dynamic approach. Ac- 
tually, the renornialization is reduced to the well known 
static renormalization group equations. In Section ^ we 
consider correlation functions of the "charge density" and 
ground the properties announced above. In Section || we 
generalize our procedure for the case of superfluid films. 
In the Conclusion we discuss the main results of our work 
and their possible relations to other systems. Some cal- 
culations are placed into Appendices. 



where the subscripts i,j label defects, positions of 

the defects, a is a cutoff parameter of the order of the 
size of the defect core, rii are integer numbers determin- 
ing the "strength" of the defects, /3 is a dimensionless T- 
dependent factor and a is the energy associated with the 
core. The expression (1.1) is correct for quantum vortices 
in superfluid fllnis, for disclinations in hexatic films, and 
for spin vortices in 2d planar magnets. For dislocations 
in crystalline films the expression (1.1) has to be slightly 
modified jsj , but the main peculiarity of the free energy, 
the logarithmic dependence on the separation, remains 
the same. 

The Gibb s di stribution exp(— JF/T) corresponding to 
the energy (1.1) can be treated as the partition func- 
tion of two-dimensional point particles with charges rij , 
(3 playing a role of the "inverse temperature" . The pa- 
rameter (3 can be considered also as the Coulomb cou- 
pling constant. Basing on the electrostatic analogy one 
can introduce the "charge density" 



p{r) 
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UjS {r — Xj) 



(1.2) 



The quantity p is vorticity for superfluid fllms and disclin- 
icity for hexatic films. We will treat the case when de- 
fects are produced by thermal fluctuations. Since both 
creation and annihilation processes conserve the "charge" 
we should accept that the total charge is zero: J^j "-j = 0. 
It leads to the constraint 



d^rp{r) = 0, 



(1.3) 



where the integration is performed over the total area of 
the specimen. 

Below we assume that for |n| > 1 the core energy fi{n) 
is so large that such defects are hardly created. Then 
only defects with the charges rii — ±1 should be taken 
into account. We will call the objects with the charges 
rii = 1 defects and the objects with the charges = — 1 
antidefects. Because of the constraint J2j '^j = there 
can be simultaneously defects and A^ antidefects in 
the system. Thus, the partition function of the system 
can be characterized via a set of probability distribution 
functions V2N depending on coordinates of 2A^ "parti- 
cles". In accordance with Eq. (1.1) the fimctions can be 
written as 



1. BASIC RELATIONS 

Static properties of the system of the vortex-like de- 
fects in thin films can be described quite universally. The 
starting point of the description is the free energy asso- 
ciated with the defects 

.F=-5]T/3n,n,ln(^l-^)+^M%), (1-1) 



V2n{xi, . . . ,X2n) 

-2-(|)™-4i:''".".'"(^)}. d") 

where Z is the sum over states and the quantity yo = 
exp(— /^/T) is usually called fugacity. The possibility to 
neglect charges with |n| > 1 implies that the fugacity is 
small. 
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The low-temperature (insulator) phase can be treated 
as a system constituted of bound defect-antidefect pairs. 
In the high-temperature (plasma) phase there are un- 
bound charges which essentially influence the system on 
scales larger than the correlation length Tc. We will treat 
the low-tempcrature phase and the region of scales be- 
tween a and rc in the high-temperature phase where one 
can neglect the role of the unbound charges and only the 
bound defect-antidefect pairs have to be taken into ac- 
count. The presence of the pairs in the system leads to 
non-trivial "dielectric" properties of the medium. As a 
result the interaction between the charges is modified, the 
effect can be described in terms of a scale-dependent "di- 
electric constant" of the medium as is suggested in Ref. 
1^. In other words, the effective coupling constant (3 be- 
comes dependent on the separation between the charges. 

The scale dependence of (3 can be described in the 
framework of the scheme proposed by Kosterlitz, [p^ . 
Namely, the partition function of the system can be in- 
tegrated over separations of the defect-antidefect pairs 
between the core size a and a scale r. After the proce- 
dure, that can be interpreted as shifting the core radius 
a ^ r, the form of the probability distribution functions 



(1.4) is reproduced (with r instead of a), but the param- 
eters P and y are renormalized. The r-dependence of 
/3 and y is determined by the following renormalization 
group equations found in Ref. jl^ 



d(3 



dln(r/a) 



-cy 



dy 



dln(r/a] 



= (2-/3)y, (1.5) 



where c is a numerical factor of order unity. The r- 
dependent function y is the renormalized fugacity. It 
determines a concentration of defects belonging to the 
bound pairs with separations of the order of r, the con- 
centration can be estimated as j//r^. The renormalized 
value of /3 determines the dependence of the strength 
of the Coulomb interaction on the separation between 
the charges. In the low-temperature phase, the effective 
value of (3 tends to a constant on large scales. The asymp- 
totic value of P is larger than 2, the critical value (3 = 2 
corresponds to the transition temperature. In the asymp- 
totic region, where (3 can be treated as r- independent, the 
renormalized fugacity y remains r-dependent. Its asymp- 
totic behavior can easily be extracted from Eq. (1.5): 



2-/3 



y (X r 



Thus, in the low-temperature phase y tends to zero as 
scale increases. 

Let us turn to simultaneous correlation functions of the 
charge density p (1.2). The odd correlation function are 



zero. Indeed, the system is symmetric under permuting 
defects and antidefects whereas the charge density ( |1.2| ) 
changes its sign at the permutation. The pair correlation 
function can be written as (see, e.g., Ref. |l9]) 



A generalization of the relation (1.7) can be obtained (see 
Ref. M) which is 



(p(ri)...p(r2„)> 



y,4n 



(p(r,)p(0))" , (1.8) 



where all separation jr^ — rj\ are assumed to be of the 
same order r*. In the large-scale limit where /3 is satu- 
rated we have 

(p(Xri) . . . p{Xr2n)) = X-2^"(p(ri) . . . p(r2„)) , (1.9) 



where X is an arbitrary factor. The relation (1.9) shows 
that the simultaneous statistics of p has normal scaling, 
that is scaling exponents of the correlation functions of 
the order 2n are equal to n times the scaling exponent 
of the pair correlation function (p_.7|). We will demon- 
strate that the behavior of non-simultaneous correlation 
functions of the charge density is quite different. 



A. Dynamics 

To examine dynamical characteristics of the system we 
should formulate a dynamical equation for a defect mo- 
tion. Following Ref. we accept the following stochas- 
tic equation 



dxj D dT 



(1.10) 



determining the trajectory of the j-th defect. Here is 
the free energy (pTl]), is a diffusion coefficient, and 
are Langevin forces with the correlation function 



(C»,a(il)0,/3(^2)) = 2DS,jSo.l3S{h - t2) . 



(1.11) 



{p{r)pm-yHr)y 



The diffusion coefficient D determines mob ility of the 
defects. We believe that the equation ( 1.10 ) is applica- 
ble to the dynamics of disclinations in hexatic films like 
membranes, freely suspended films and Langmuir films. 
The equation for the vortices in superfluid films is a bit 
more complicated. It is written in Section ^ where the 
correlation functions of the vorticity are analyzed. 

The equations ( 1.10| , m| ) describe trajectories of sep- 
arate defects. We shou d also take into account annihila- 
tion and creation processes. Remember that we neglect 
(1.6) defects with \nj\ > 1. Next, processes where a number of 
defect-antidefect pairs are created at the same point are 
suppressed since probability of such events is small due to 
the energy associated with the cores of defects. Then we 
have to take into account the creation processes of single 
pairs solely, they are characterized by the creation rate 
R{r) which is a probability density for a defect-antidefect 
pair with the separation r to be created per unit time per 
unit area. The annihilation processes have to be charac- 
terized by the annihilation rate R{r) which is a proba- 
bility for a defect-antidefect pair to annihilate per unit 
time if the pair is separated by the distance r. Really, 



(1.7) 
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both R{r) and R{r) are nonzero only if r is of the order 
of the core size a. Let us introduce the integrals 



(1.12) 



A= / d^rR{r) 



A = / dVi?(r) 



Here, the creation constant A is a probability for a defect- 
antidefect pair to be created per unit time per unit area 
and A is a constant having the same dimensionality as the 
diffusion coefficient D. Below, the diffusion coefficient D 
is put to unity by rescaling time. Then the annihilation 
constant A is a dimensionless parameter of order unity 
and the creation constant A can be estimated as 

A-a^''exp(-2^/T), (1.13) 

which is the second power of the defect concentration. 

d'P2N _ 

at 



As in statics, in 
can be described in 
functions. In our 

the system is zero, 

V2Nit, Xi,. . . , Xn, Zi, 

and Zj are positions 
fects correspondingly, 
equal to unity which ^ 



dynamics the system of defects 
terms of probability distribution 
case, when the total charge of 

only even probability densities 
...,zn) are non-zero, where Xj 
of the defects and of the antide- 

The total probability should be 
nves the normalization condition 



1 

xV2N{t,Xi, 



d^xi . . . d'^xjy d^zi . . . d^z^r 



, zn) 



1. 



(1.14) 



Starting from the equation ( 1.10 ) and taking into account 
the creation and annihilation processes one can derive a 
system of master equations for the probability densities 





' 1 ^-^ ^ ' 


d 


'1 dT^ ' 


' 3 




[dxj 




dzj 




dx'j dz'j 



^R{xj - Zk)V2N + J d'^xd'^zR{x - Z)P2N+2{XI, . . . ,Xn,X,Zi, . . . ,zn,z) 

Xn, Zi, . . . , Zk-1, Zk+1, . . . , Zn) — XA'P2N , 



(1.15) 



^R{Xj - Zk)'P2N-2{xi 



, Xj — i , Xj-^i 



where A is the area of the film. 



The Gibbs distri butio n (1.4) must be a solution of the 
master equations (1.15). The condition imposes the fol- 
lowing constraint on the creation and the annihilation 
rates 



i?(r) 



(1.16) 



where we imply r > a. The constraint ( 1.16 ) can be 
treated as the manifestation of the equilibrium state of 
the thermal bath which in our case is related to short- 
scale fluctuations. Thus the constraint ( 1.16| ) has the 
same origin as Eq. ( 1.11 ). At deriving Eq. ( 1.16 ) we 
assumed that the separation |a; — 2;| in the argument of 
the rate i? or ^ is smaller than separations between x 
or z and other points. That is accounted for the small 
characteristic value of the separation which are of the 
order of the core radius a. N ote t hat for such r the fac- 
tor yQ{r / af^^^ entering Eq. (1.16) can be treated as the 



r cno rmalized value y of the fugacity as follows from Eq. 

(ph. 



Principally, the master equations ( 1.15 ) enable one to 
find conditional probability densities, related to different 
time moments. Consequently, starting from the equa- 
tions one can examine non-simultaneous correlations in 
the system. However, there are terms in the master equa- 
tions (associated with the creation and annihilation pro- 
cesses) mixing the probability densities V2N with differ- 
ent N . That makes the master equations hardly useful, 
that is a motivation to look for some more suitable tech- 



nique. Such a technique was developed by Doi, Ref. 
we formulate it in the subsequent subsection. 



B. Quantum Field Formulation 

The Doi technique ||l^ enables one to treat systems of 
classical particles where creation and annihilation pro- 
cesses occur. The main idea introduced by Doi is that 
correlation functions of different quantities characteriz- 
ing the particles can be written in the form close to the 
one known in the quantum field theory. Of course there 
are some peculiarities related to the fact that for classi- 
cal particles one should deal directly with probabilities 
whereas in the quantum field theory one starts from the 
scattering matrix. Nevertheless the Doi technique en- 
ables, say, to formulate a diagrammatic expansion with 
the conventional rules. The technique was originally de- 
veloped to describe systems of molecules involved into 
chemical reactions. But it is definitely applicable also to 
the system of point defects. 

The Doi technique is formulated in terms of the cre- 
ation ■(/; and annihilation ip operators which satisfy the 
same commutation rules as ones for Bose-particles 



[il){ri),i){r2)] = <5(ri - ra) , 



(1.17) 



For our system of defects we should introduce annihila- 
tion and creation operators '0± and ^± where the sub- 
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scripts + and — label fields related to the defects and 
to the antidefects. The state of the system at a time 
moment t can be written in terms of a "quantum" state 



1 



N=0 



(N 



1^2 



d xi . . .d xn d zi . . . d zj^V; 



2N 



x-0+(a;i) . . . ■4!^{xn)'iP-{zi) . . . -0_(2:7v)|O) 



(1.18) 



where V2N are the probability densities introduced above 
and |0) designates the vacuum stat e: V'ilO) = 0. In ac- 
cordance with the expression ( |l.l^ ) an evolution of the 
quantum state \t) is determined by the master equations. 
The evolution equation can be written as 

dt\t)^-n\t), and |t2> = exp [(ti - t2)W] , (1.19) 

where TL is an operator expressed in terms of the fields 
V^it and '0± • By analogy with the quantum field formula- 
tion it can be called the Hamiltonian operator or simply 
the Hamiltonian. 

The total probability must be equal to unity, which 
leads to Eq. (1.14). The condition can be written in 
terms of the quantum state \t) as 



(sum|t) = 1 , 



(1.20) 



where 

(sum| — (0| exp 
Note the following identity 



dV (-0+ +-0-) 



(sum|0±(r) — (sum| , 



(1.21) 



which can be easily checked using the commutation rules 
( 1.17 ) and the equality {Q\ip± = 0. Since the evolution 
must conserve the total probability (sum|t) the relations 

(sum|7i = , and (sum| exp {tH.) ~ (sum| , (1-22) 

have to be satisfied, where r is an arbitrary parameter. 

Quantities characterizing the system can be repre- 
sented by corresponding operators, see Ref. [^. Say, 
the operator of the charge density is 



(1.23) 



If A is such operator corresponding to a quantity A, then 
an average value of the quantity at a time moment t can 
be expressed as 



Note that 



{A{t)) = (sum|i|t) . 
(■0±) = (sum|?/>±|t) = 1, 



(1.24) 



(1.25) 



whic h is a consequence of Eqs. (1.20,1.21). The rela- 
tion (1.25) shows that it is natural to shift the creation 
operators ip± introducing new variables (see Ref. pl||) 



0± = 1 + 0± : 



0. 



(1.26) 



Correlation functions of dif feren t quantities can be pre- 
sented analogously to Eq. ( 1.24 ). For example, the pair 
correlation function of two quantities ^(^i) and B{t2) 
(we assume ^2 > ^i) can be written as 

{B{t2)A{h)) = (suni|Sexp [(^i - t2)H] A\ti) . (1.27) 



Usin g the relations ( |l.l9[|l.22| ) we can rewrite the expres- 
sion (1.27) as 



{B{t2)A{h)) - (sum|S(i2)i(ii)|in), 



(1.28) 



where A{t)^B{t) in Eq. ( |1.28D are operators in the 
Heisenberg representation 



A{t) = exp [~{ti - t)H] icxp [-{t - tin)n] 



(1.29) 



satisfying the equation dtA{t) = A{t) . In Eq. ( |l.28| ) 

I in) is an initial state (realized at a time moment tin) and 
tf is a "final" time, so t hat tf > t2 > ti > tin- 

The expressions like ( 1.28| ) enable one to reformulate 
the problem of calculating correlation functions in terms 
of a functional integral, see Ref. [Q. Namely, we can 
write 



{Aiih) . . . An{tn)) 



X exp 



Vij±Vtp±Ai ...A„ 



+ / dV [V'+(tf,r) + 0_(tf,r)] L 



(1.30) 



where ^± , 'ip± are to be interpreted as functions of t and 
r. We assume that tf > ti,. . . ,tn in Eq. ( 1.30 ). Deriv- 
ing the expression one has taken the limit tin — *■ — c» and 
assumed |in) = |0). Because of the creation processes the 
vacuum has to be turned into a stationary state during 
the infinite time. To ensure convergence of the functional 
integral (1.30) the integration contour over the field -ip 
shoul d go parallel to the imaginary axis. Note that the 
shift ( 1.26 ) kills the boundary term J d^r {ip++ip-): It is 
cancelled by a contribution originating from the deriva- 
tives 9tV'± after integrating over time. Then we come to 
a conventional representation of the correlation functions 
in terms of a functional integral 



(Ai(ti)...A„(t„)) 



'D%p±'D'ip± exp 



dt 



Ai...Ar, 



n 

(1.31) 



The relation (1.30) or (1.31) is a convenient starting point 
for treating a system of particles involved into the cre- 
ation and annihilation processes. 
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2. DIAGRAMMATIC REPRESENTATION 

Below, we apply the Doi technique to our particular 
problem. The explicit expression for the Hamiltonian 
determining the e volut ion of the defect system in accor- 
dance with Eg. ( 1.19| ) can be found from the master 
equations (1.15). Comparing the equations with Eqs. 
( priSplgj ) we get 



P ■ 



(2.1) 



The explicit expressions for the terms entering Eq. (2.1) 
are 



Hr = - / dVidVa 



d^r V?/;+V-0+ + V-0_Vi/;_ 



R{ri - r2)(V'+,i^-,2 - 1) 



(2.2) 



+i?(ri - r2)(V'+.l'0-,2 - '0+,l'0-,2V'+4V'-,2) 



(2.3) 



ri - r2 



\ri - r2 



-?/;+4V-,2 



2/3 / dVi dVa 



\ri - r2 



(2.4) 



where V'+.i — '0+ (^,ri) and so further. The diffusive 



contr ibu tion (2.2) is related to the Langevin forces, in 
Eq. (2^) -R is the annihilation rate and R is the creation 
rate for the defect- antidefect pairs (th e qu antities were 
introduced in Section]^), and the term ( |2.4|) des cribe the 
Coulomb interaction. Using t he pr operty (1.21), one can 
easily ch eck the conditions ( 1.22|) for all contributions 



(PP) 



Performing the substitution ( 1.26| ) we can express the 
Hamil tonia n (2.1) in terms of the fields -i/Ji. Note that the 



shift (1.26) kills terms of the second order proportional 
to A, A and generates additional third-order vertices. Of 
course one can work in both representations. It is more 
convenient for us to use Eqs. (1.30,2.1). We can eas- 



ily convince oursel ves th at odd correlation functions of 
the char ge de nsity ( 1.23 ) are zero. Indeed, the exponent 
in Eq. (1.30) is invariant under permuting ^ 

<-> ip-j whereas the charge dens ity changes its sign 
at the permutation. The constraint (1.3) shows that the 
symmetry is not spontaneously broken what could lead 
to non-zero odd correlation functions. 
It follows from Eq. 



1.29) that the commutator [H,p\ 



should be equal to — Vj where j is the current density op- 
erator. Calculating the commutator with Eqs. (1.23,^) 
we get 



-(V3+^++V^-V^-)V0, (2.5) 



where (f> is an "electrostatic potential" 



(r) = -2/3 / d'x In 



\r — x\ 



~p{x). (2.6) 



Principally, besides the "internal" potential ( |2.6| ) an "ex- 
ternal" potential (/)ext can be imposed onto the system, 
satisfying the equation V^^cxt = 0. Then an "exter- 
nal force" should be added to the right-hand side of the 
equation (1.10). The force generates an "external" con- 



tribution to the Hamiltonian 



He 



d^r 



(2.7) 



The expression (2^) for the current density has also to 
be corrected by substituting (j) ^ (j) + (j>e^t- With the 
term (2.1) we can examine susceptibilities describing a 
response of the system to the external influence. 



Substituting the expression (2.1) into (1.30) or (1.31) 
and expanding the exponent over Tin and Tip one can 
obtain a conventional perturbation series for calculating 
different correlation functions of Tp,ip- The series is an 
expansion over R, R and (3 in terms of the conventional 
diffusion propagators: 

G(i,r)-(7/i+(i,r)^+(0,0))o 

= (^_(t,r-)^-(0,0))o = l^exp ("-Q , (2.8) 



where 6{t) is the step function. However, effects related 
to the Coulomb interaction and to the annihilation pro- 
cesses are not weak. Therefore one must take into ac- 
count the Coulomb interaction and the annihilation pro- 
cesses exactly. By other words, at calculating the cor- 
relation functions one must consider the complete series 
over P and R. Fortunately, the the expansion over R is 
equivalent to an expansion over the fugacity y which is 
assumed to be a small parameter. Therefore we can take 
only principal terms in the expansion over R. 

The perturbation expansion can be formulated as a di- 
agrammatic series. We develop the diagrammatic tech- 
nique starting from the representation (1.3C), pushing 
the final time tf to the far future. We depict the prop- 
agator (2^) by a line directed fr om ip to ip. The term 
with the creation rate R in Eq. (^.l[) generates vertices 



where two propagator lines start, the vertices correspond 
to the defect- antidefect creation processes. The Coulomb 
term in Eq. (|]^) generates two-point objects which we 
will designate by dashed lines, such line describes the 
Coulomb interaction of defects located in points con- 
nected by the line. And th e te rm proportional to the 
annihilation rate R in Eq. (2.1) produces two types of 
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vertices. First, it produces vortices where two propa- 
gator lines finish, that corresponds to an annihilation 
process. Second, it produces fourth-order vertices which 
correspond to an effective interaction related to a finite 
probability for a defect-antidefect pair to annihilate, see 
Ref. jl^ . A typical diagram block is presented in Fig. ||. 
The block is drawn in real r — t space-time. The curves 
constituted of the propagator lines can be interpreted as 
trajectories of defects and antidefects. Due to causality 
the particles always move forward in time. Note that the 
dashed lines corresponding to the Coulomb interaction 
are perpendicular to the t-axis since the interaction is 
simultaneous. 



A. Pair Conditional Probability 

In the subsection we examine an auxiliary object which 
will be needed for us at intermediate stages of subsequent 
calculations. The object is the following correlation func- 
tion 

M{t2 - ti,ri,r2,r3,r4) 

= (V'+(t2,r-i)V'-(t2,r-2)^+(ii,r-3)^-(ii,r4)) . (2.9) 

For a stationary case the average ( |2.9| ) depends on the 
difference t — t2 ~ ti only. Due to caus ality M is equal 
to zero provided t < 0. The quantity ( ^.9[ ) can be in- 
terpreted as a probability density to find a defect and 
an antidefect at the time moment t2 in the points ri 
and r2 provided they were located in the points ra and 
at the time moment ii. It can be considered also as 
a two-particle matrix element of the evolution operator 
exp[-(t2 - ti)TC\. 

As we explained above, the perturbation series in terms 
of the creation rate R is an expansion over a small param- 
eter. Here we examine the principal contribution to the 
conditional probability (|2.9|) w hich is of the zero order 
over R. Then the average ( |2.9| ) can be represented as a 
series of diagrams of the type depicted in Fig. ^. One can 
interpret the picture as trajectories of a defect and of an 
antidefect which are driven by the Langevin forces, and 
are influenced the Coulomb interaction (dashed lines) 
and the effective interaction associated with the anni- 
hilation processes (point vertex). Note that in this ap- 
proximation direct annihilation events do not contribute 
to the conditional probability (2.£) since they would lead 
to terminating the lines in the diagrams. 

It is of crucial importance that both the Coulomb in- 
teraction and the effective interaction associated with the 
annihilation processes are local in time. Therefore all the 
diagrams representing the conditional probability ( |2.9| ) 
are ladder diagrams, like in Fig. ^. Summing up the 
ladder sequence we get an equation for M which can be 
written in the differential form 



dtM = (V? + V^) M + 2/3 (Vi - V2) 



M 



.2 , „2^^^^2/3(Vl-V2) ^^^^ 
-i?(ri - r2)M + 5{t)S (ri - rg) 5 (r2 - r^,) . (2.10) 



Since M = at t < we conclude from Eq. ( ^.10 ) that 
at i -1-0 

M{t, ri,r2, ra, r4) -> ^ (ri - r^) S (r2 - r4) . (2.11) 

The total probability to find the defect-antidefect pair 
is determined by the integral / d^ri d^r2 M. Let us cal- 
culate the time derivative of this integ ral substituting 
dtM by the right-hand side of Eq. ( |2.iq ). Then the first 
two terms will give zero contributions (since they are to- 
tal derivatives) and only the term with R will produce 
a non-zero (negative) contribution. It is quite natural 
since the Langevin forces and the Coulomb interaction 
cannot change the total probability whereas the annihi- 
lation processes diminish it. Note that all the terms in 



the right-hand side of Eq. ( 2.10 ) proportional to M have 
the same dimensionality. Therefore one could expect a 
simple scaling behavior when t scales as r^. The subse- 
quent calculations confirm the expectation. 
In terms of the variables 



ri - ^2 , 



^"1 + f2 



r-o = r3-r4, (2.12) 



the equation ( 2.10 ) for M is rewritten as 
'1 



dtM 



-yt + 2V,^ + 



R{r) M 



+6{t)S (r - ro) S (g - r^j^ - r4/2) . 

We see that the differential operator in the right-hand 
side of the equation falls into two parts depending on 
Q and r only and that the "source" is a product of 5- 
functions of the same variables. Therefore the solution 
of the equation can be written in a multiplicative form 



the function S satisfies the following equation 



dtS^^^'-S 



/ r 



S 



-R{r)S + 6{t)5{r-ro) 



Note that 



if then S ^ 6 {r - Vq 



(2.14) 



(2.15) 



The relation follows from Eq. (2.14) and causality (lead- 
ing to S' = for negative t). 

In accordance with Eq. ( ^.13 ), a motion of the mass 
center and the relative motion of the defects are sepa- 
rated. The motion of the mass center is purely diffu- 
sive whereas the relative motion is strongly influenced 
by the interaction. The function S can be treated as the 
probability density for the relative motion of the defect- 
antidefect pair. It is natural to expand the function into 
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the Fourier series over the angle between the vectors r 
and Tq: 



+00 



S{t, r, To) = ^ Sm{t, r, ro) ex.p{im(p) . 



(2.16) 



Motions corresponding to different angular harmonics are 



separated. In terms of the angular harmonics Eq. (2.14) 
is rewritten as 



1 



1 j-f-) ^ 

a2 + (i + 2/3)-a„- — 



■<5(t)<5(r - ro) 



4:7Tr{ 



(2.17) 



I t is possi ble to get equations for S analogous to Eqs. 
(2.14 2.17) in terms of rp. They have practically the 



same form as Eqs. (2.14,2.17). The only difference is in 
the sign of (3 which is opposite. That leads to the relation 

Sra{t,r,ro)= {^^y S^it,ro,r). (2.18) 



Let us stress that the relation ( 2.1g ) is correct for an 
arbitrary function R{r). 

Consider a behavior of the angular harmonics 
r, ro) at small r. More precisely, we assume f> 
and examine the regi on V 7 ^ r ^ a. Then it is possible 
to use the equation ( |2.17| ) with the time derivative and 
the annihilation term neglected. As a result we get 



2v 



(2.19) 
(2.20) 



where Ci_m,C2,m are some factors dependent on t and 
ro . The ratio of the factors is determined by a concrete 
r-dependence of the annihilation rate i?, one can assert 
only that Cim and C^m are of the same order. Therefore, 
if we consider the behavior of the function S for r ^ a , 
then the second term in the right-hand side of Eq. ( 2.19| ) 
can be neglected. By other words, being inter ested in the 
scales r a, we can solve the equation ( ^.17| ) neglecting 
the annihilation term and requiring a finite value of Sm 
at r ^ instead. 

The requirement can be treated as the boundary con- 
dition for Sm at small r. The other boundary condition 
is that Sm tend to zero at r — > 00. The equations fo r 5*,^ 
with the boundary conditions are solved in Appendix A 1 . 
We present here only the answer 



Sm{t, r,ro) 



1 /roy 



exp 



8t 



\ At ) 



(2.21) 



where / is the modified Bessel function and 1^ is intro- 



duced by Eq. (2.20). Remind that the expression is cor- 
rect provided r, ro,-\/t 3> a. Extracting from Eq. (2.21) 
an asymptotics at small r we get 



Cl.m — 



87rr(l + iy)t 



exp 



10. 

St. 



Note that the Coulomb term in Eq. ( 2.14 ) produces 
a probability flux to th e orig in. To find it we should 
integrate the equation ( 2.14 ) over a disk of a radius 
a <^ r <^ ^/i centered at the origin and single out the 
contribution to dt J (PrS associated with the Coulomn 
term. Then we find the flux XrC'1,0 where 



Xr = 8n(3 . 



(2.22) 



One can treat the quantity (2.22) as the renormalized 
("dressed") value of the annihilation constant. Now we 



understand why the solution (2.21) (realized at r ^ o) is 
insensitive to a particular form of the annihilation rate. 
The probability for a defect-antidefect pair with the sepa- 
ration r 3> a to annihilate is determined by the Coulomb 
attraction. And only the behavior of the probability den- 
sity at r ^ a is sensitive to the particular form of t he an - 
nihilation rate i?(r): The coefficients C2,m in Eq. ( 2.1£ ) 
are positive if A < Ar and are negative if A > Ar . 

Returning to the conditiona l probability we ob- 
tain from Eqs. (|2.13|,^.16p^) 



M = 



X exp 




exp(im(^) I, (^) , (2.23) 



m— — 00 



\ At J 



where v is introduced by Eg. (2.20). One can easily 
check that the expression ( 2.23 ) is reduced to Eq. ( 2.11 ) 
at i — > -1-0. It is possible to calculate an explicit ex- 
pression for the total probability to find a defect and an 
antidefect in a ny points at a fixed time separation t, see 
Appendix A 1 , The asymptotic behavior of the expres- 



sion (2.23) at the condit ion r rp/t 3> 1 is examined in 
Appendix [A 2| , see Eq. (A14). If both r, ro are much 
greater than \/t, it can be written as 



(2.24) 




X exp 



The answer is quite natural. The characteristic values of 
the separations ri — and of r2 — r4 are of the order 
of y/t and are consequently much smaller than |ri — 
(or |r3 — r4|). Then r « ro and it is possible to neglect 
all terms, containing |ri — in the denominators, in 
the equation (2.10). Thus we come to a pure ly diffusive 
equation leading to the asymptotic law ( 2.24 ). 

The consideration presented above can be generalized 
for the case when an "external electrostatic potential" 
(j)cxt is imposed onto the system. Its infiuence is described 
by the contribution (2.7) to the Hamiltonian. Performing 
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the same procedure as abov e we get a modified equation 
for the correlation function (2.9) 



Vl) M 



4f3\/r -^M - R{r)M 



+V^ext(t2,ri)ViM- V(/)cxt(t2,r-2)V2M 

+5{t2 - h)6 (ri - rg) S (rs - Vt) . (2.25) 

The expression ( 2.25| ) shows that the gradient of the ex- 
ternal potential has to be added to the gradient of the 
internal one. Of course, in the pres ence of the exter- 
nal field the correlation function (2.9) depends on both 
time moments ti and t2. Note that the operator in the 



right-hand side of Eq. (2.25) is the same as that for the 
Fokker-Plank equation formulated in Ref. [0 (excluding 
for the annihilation term). 



3. RENORMALIZATION 

In this section we are going to discuss effects related to 
high-order terms over the creation rate R. The effects are 
relevant only near the transition point where (3 is close 
to 2. Then the influence of small-scale defect-antidefect 
pairs on larger scales becomes essential. In the situation 
the most natural language is the renormalization group 
approach. One can formulate a renormalization group 
procedure in the spirit of Kosterlitz, Ref. We will 

single out blocks corresponding to small separations of 
the pairs and treat the m as renormalized quantities en- 
tering the Hamiltonian (2.1). 



A. Creation and annihilation rates 

Sizes of the pairs are small near creation and near anni- 
hilation points. Here, we consider vicinities of the points. 
Then it is possible to neglect the interaction of the de- 
fect and of the antidefect with the environment. Thus 
we turn to the situation when only a single pair can be 
treated. If this is the case then one should analyze dia- 
gram blocks of the type drawn in Fig. |. The left part of 
the figure corresponds to a vicinity of the creation occur- 
ring at a time moment ti and the right part of the figure 
corresponds to a vicinity of the annihilation occurring at 
a time moment t^. 

Consider processes occurring during a time interval 
T from the creation time ti. One can separately treat 
a block corresponding to the time interval from ti till 
t2 = ti + T. For the purpose we use the well-known 
property of the propagators ( |2.8| ) 

G{s3 - si,r) = J (fxG{s3 - S2,r - x) 

xG(s2-5i,a;), (3.1) 

where S3 > S2 > si. For each diagram we extract propa- 
gators G containing t2 inside their time interval and rep- 
resent the propagators like in Eq. (3.1) beheving S2 — t2- 



The procedure is reflected in Fig. || where the dotted line 
represents a plane t = t2 ui the r — t space-time and the 
integration in Eq. (3.1) corresponds to the integration in 
the plane. As a result, the block to the left of the plane 
is separated, it is characterized by the time separation 
T and by two points ri and r2 lying in the plane, the 
points are intersections of the plane with the trajectories 
of the particles. The block has to be inserted into more 
complicated objects via a convolution over ri and r2. 

The same is true for the vicinity of the annihilation 
point also. Let us take a time moment separated by a 
time interval r 3> from an annihilation time t^. Then 
it is possible to introduce the block which is a sum of 
the diagrams where the trajectories of the annihilating 
particles start from two given points and at t — t^. 
The block has to be inserted into more complicated ob- 
jects via a convolution over the points. In the vicinity of 
the annihilation point we can take into account the in- 
teraction of the annihilating defect-antidefect pair solely. 
That leads to the same ladder diagrams treated in Sec- 
tion ^ Therefore we can write an expression for the block 
without an additional analysis 

Rr{ro)= J dVi dV2i?(r)M(T,ri,r2,r3,r4) 

= J d^rS{T,r,ro)Rir). (3.2) 



Here r and Tq are defin ed by Eq. ( ^.12 ), M is the condi- 
tional probability (2.9), S is the conditional probability 
for the relative motion of the defects, see Eq. (2.13), it is 
the solution of the equation (|2.14). The physical meaning 



of the quantity Rriro) is a distribution of the annihilat- 
ing particles over the separation rg between the particles 
at the time moment t^. It is natural to name this dis- 
tribution "dressed" annihilation rate since the quantity 
determines a probability for the particles to annihilate 
after the time interval r. Note that all processes occur- 
ring on scales larger than ^/T are sensitive only to this 
dressed quantity. 

Let us substitu te in to Eq. (3.2) the product RS ex- 



pressed from Eq. ( p.l4| ) . The terms with the total deriva- 
tives give zero contribution to the integral over r and we 
get 



i?r(ro) = -dr / d'rS{T,r,ro), 



(3.3) 



Since S{t) tends to zer o at r - 
negative r we get from ( ^.15 3.3 ) 



J drRriro) = 1. 



and is zero for 



(3.4) 



The relation means that the total probability for a given 
pair t o annihilate is equal to unity. As is seen from Eq. 
( ^.19 ) at the condition t ^ the main contribution to 
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the integral in the right-hand side of Eq. (3.2) is asso- 
ciated with the region r ~ ^Jt and therefore the contri- 
bution to the integral associated with the region r ^ a 



is negligibl e. Th erefore we can use the expressi on ( 2.1i ) 
with Eg. (2.21). Substituting the expression (A3) into 
Eq. (3^) we get a universal expression for the dressed 
quantity Rriro) which is insensitive to the bare quantity 

m- ^ ^ 

In Sec. we established the renormalized value ( 2.22| ) 
of the annihilation constant A. This analysis concerned 
the fourth-order interaction term written in Eq. (2.1). 
Below we demonstrate that the renormalized coefficient 
at the second-order annihilation term has the same value, 
i ndep endent of the bare one. In accordance with Eq. 
( 1.12 ), to find the renormalized value Ar we should cal- 
culate the i ntegral of Rr{ro). As is demonstrated in Ap- 
pendix Al_ at r ^ the value of the integral is inde- 
pendent of T and coincides with the value written in Eq. 
22), as one anticipated: 



(1 



K = / d^roi?.(ro) =87r/3. 



(3.5) 



The phenomenon resembles the renormalization of the 
reaction rate due to diffusion, see Refs. p^j2^ . 

Analogously, one can introduce the renormalized cre- 
ation rate Rr{r) which is determined by the block de- 
scribing the vicinity of the creation point (see Fig. ||). 
Summing up the same ladder sequence of the diagrams 
we get 



Rr{r)^ / d2r3d2r4i?(ro)M(T,ri,r2,r3,r4) 



d2roi?(ro)S'(T,r,ro) . 



(3.6) 



Here r and Tq are defined by Eq. ( ^.12 ), M is the condi- 
tional probability (2.9), S is the conditional prob abilit y 
for the relative mot ion of th e defects, see Eq. _( 2.13 ). 
Using the relations ( Ll£ , 2.18|) we get from Eq. (^^6^ 



i?.(r)=4(-)^^^«. 



(3.7) 



Thus we see that the relation (1.16) is reproduced for the 
renormalized quantities i?r and Rt. 

The renormalized creation rate Rt{t) can be inter- 
preted as a probability density to find a defect-antidefect 
pair with a space separation r provided the pair was 
born on time separated by r from the measurement. Let 
us calculate the total probability density Ar to find the 
defect-antidefect pair at a fixed time separation t regard- 
ing T ^ a^. The probability is determined by the inte- 
gral of R rir) over r. We conclude from the expressions 
(2.21 3^) that the integral is determined by the region 
r ~ ^/t. Taking into account Eq. (|3.5|) we get 



d^rRrir) 



(3. 



We see that due to annihilation of defects at collisions the 
total probability diminishes at increasing the time sepa- 
ration T as a power of r. The property can be interpreted 
as follows; The majority of defect-antidefect pairs anni- 
hilate fast after their creation and only a minor part of 
the defects achieve a separation r ^ a. The probability 
of such event is proportional to [r/a)^'^^. 

The results obtained in the subsection are correct if the 
variation of the coupling constant j3 on the scale interval 
a < r < ^/t is small. The existence of such interval is 
justified by the assumed small value of the fugacity j/o- 
Near Tc variations of /3 on a wide region of scales can be 
relevant. Then the consideration needs a generalization 
made in the last subsection of this section. 



B. Coulomb interaction and diffusion coefficient 

Let us consider the renormalization of the Coulomb 
interaction related to small defect-antidefect pairs. It is 
known that the influence of such pairs can be described 
in terms of a contribution to the effective dielectric con- 
stant, see Ref. The picture is naturally generalized 
for the dynamics. 

Before proceeding to calculations, it wil l be convenient 
for us t o ex press the Coulomb part (2.4) of the Hamil- 
tonian ( ^.l[ ) in an alternative form. Namely, using the 
Hubbard-Stratonovich trick we rewrite the fourth-order 
term T-Lp as a functional integral over auxiliary fields a 
and (b 



exp — / AtTLp 



'D(j) Va exp 



dt {Hi + n2 



Hi = / d^r 



d^r VctV0 



^2 = 



47r/3 



(3.9) 

(3.10) 
(3.11) 



The relation ( |3.9| ) can be easily checked using the bare 
expression 



\ri - r2|^ 



(3.12) 



following from Eq. (3.11). We see that the correlation 
function ( |3.12 ) corresponds to the dashed line on the di- 
agrams. Note that the field (f> in the expr essi ons is the 
electrostatic potential introduced b y E q. (2.6). Indeed, 
integrating over the field a in Eq. (3.9) we get the Pois- 
son equation 



P, 



(3.13) 
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leading to the expression (2.6). 

Now we can work in terms of the sum Tio + + 
7^1+2^2 where two first terms are introduced by Eqs. 
(2.2 2.2). Note that the sum Hq + Hi is invariant under 



the foUowing infinitesimal transformation 

(SV'h- = aip+ Sip-^ — —aip- Sip+ — —aTp-\ 
Sip- — a?/'- S(l> — 2a 5<j — V^a , 



(3.14) 



where a is a function of coordinates. If to substitute 
R{r) X6{r) and R{r) X5{r) into the expression 



(2.3) then it will be invariant under the transformation 
( 3.14 ) also. Deviations from the symmetry related to r- 
depcndencics of R and R are irrelevant. If the function 
a contains only terms linear and quadratic over r then 
the contrib ution Ti.2 (3.11) is invariant under the trans- 
formation (3.14) also. The symmetry leads to a number 
of the Ward identities. Particularly, they connect the 
renormalized triple vertices to the self-energy function of 
the propagator. 

A typical diagram contributing to renormalization of 
the effective "dielectric constant" is drawn in Fig. ||. 
There we see a loop composed of the trajectories of a 
defect and of an antidefect which annihilate after their 
creation. There are also two "external" dashed lines cor- 
responding to the interaction of the defect-antidefect pair 
with an environment. Besides the diagrams of the type 
drawn in Fig. |^ there are also diagrams with two external 
dashed lines attached to the same trajectory. We draw 
the external lines with arrows to remember that two sides 
of the dashed line representing the correlation function 
( |3.12 ) are not equivalent. We imply that the dashed lines 
are directed from the field a to the field (p. 

As previously, we can dissect the diagram into parts 
which can be treated separately. Then the answer can 
be found as a convolution of the corresponding expres- 
sions. We perform the dissection along the planes in the 
r — t space-time perpendicular to the i-axis and corre- 
sponding to the time moments t2 and tj, of the external 
Coulomb lines. In Fig. || the dissection is shown by the 
dotted lines. We see that the loop is divided into three 
parts. 

The left part of the loop implying the integration over 
the time ti (see Fig. ||) corresponds to 



dr / d rsd r4 i?(r3 - r4)M(r, cci, a;2, ^3, r4) 



AtRt{\xi - X2\ 



(3.15) 



Substituting here Eq. (|3.7D and using Eq. (3.4) we get 



2/3 



(3.16) 



The central part of the diagram depicted in Fig. g 
corresponds to the conditional probability (2.9) M(<3 — 
t2, a;3, a;4, a;i, 0:2). And the right part of the diagram in 



Fig. ^ corresponds to the integral (3.4). The relation 
can be recognized as a manifestation of independence of 
all results of the final time in the relation (1.30). If 
we chose if = then the right part of the diagram in 
Fig. H disappears and we should substitute 1 instead, in 
accordance with (3.4). 

The structure of the diagram depicted in Fig. ^ shows 
that the block related to the defect-antidefect pair can be 
treated as a self-energy insertion to the line correspond- 
ing to the Coulomb interaction. Thus it is natural to 
expect that this insertion can be treated as a contribu- 
tion to the "dielectric constant" , leading to a renormal- 
ization of the Coulomb constant (3. The corresponding 
quantitative analysis is presented in Appendix ^ giving 
the following expression for the correction to /3 



A/3: 



a\2/3-4 



(3.17) 



The expression can be treated as an integral over the 
characteristic sizes of the defect-antidefect pairs. 

One may try to find more complicated blocks con- 
tributing to a renormalization of the Coulomb coupling 
constant (3. An example of such block is depicted in 
Fig. ^ where a number (three) "external" lines are at- 
tached to the loop corresponding to the trajectories of 
the defect-antidefect pair. One can easily check that the 
block depicted in Fig. gives a correction to the Coulomb 
force which diminishes faster than at increasing the 
distance r between the interacting particles. Therefore 
the contribution is irrelevant. The same is true for more 
complicated diagrams of the same type. 

We can also consider blocks which can be treated as 
contri butions to the diffusion coefficient D introduced by 
Eqs. (1.10,1.11). An example is depicted in Fig. ^, where 
the block between two dotted lines is a self-energy inser- 
tion to the propagator (2.8) which gives the renormaliza- 
tion of the diffusion coefficient. We will assume that the 



fields V'± are corrected to keep the term (2.2) unchanged. 
Then the contribution ( 3.18 ) has to be extracted from 
the renormalization of the coefficient in front of the time 
derivatives. 

To analyze the correction AD quantitatively one 
should know a three-particle conditional probability 
which is more complicated than the two-particle condi- 
tional probability ( |2.9| ). Fortunately, one can estimate 
the value of AD without detailed calculations. The point 
is that the dependence of AD on the cutoff a can be pro- 
duced only by regions near the creation or near the anni- 
hilation point (which are designated by ovals in Fig. ||). 
The regions can be anal yzed in terms of the two-particle 
conditional probability (2.9) since only the interaction 
of the nearest "particles" is relevant there. We already 
know the answer: the region near the creation point 
produces the renormalized creation rate ( |3.6D whereas 
the region near the annihilation point produce no a- 
dependence. Then simple dimensional estimates give the 
answer similar to the expression (3.17) 
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AD 



-yl 



dr /a\ 2/3-4 
r 



(3.18) 



Remind that the bare value of D is assumed to be equal 
to unity. 

Of course there are also blocks which can be inter- 
preted as corrections to the triple vertices describing the 
interaction of the fields tp± , tjj± with the Coulomb fields 
(j) and a. An example of such block can be imagined if 
to attach an external dashed line to a trajectory between 
the dotted lines in Fig. ^. However, in the renormaliza- 
tion scheme accepted the corrections are fully absorbed 
into the renormalization of the fields ip± . That is a con- 
sequence of the symmetry of the Hamiltonian under the 
transformation (3.14). Namel y, th e cubic term ( 3.10D is 
unchanged provided the term (2.2) is unchanged. 



C. Summary 

In the previous su bsection w e obtained the expressions 
for the corrections (3.17,3.18) of the Coulomb constant 
and of the diffusion coefhcient in the main order over 
the fugacity t/o- Now we are going to discuss high-order 
corrections over uq whic h in statics lead to the renormal- 
ization group equations (1.5). It will be more convenient 
for us to proceed in spirit of the Kosterlitz renormal- 
ization grou£_scheme. Namely, we see that the expres- 



sions ( 3.17 , 3.181) are written as integrals over the space 
variable r which can be treated as the size of a defect- 
antidefect pair. We can first perform the integration over 
a restricted interval of the sizes, what gives slightly renor- 
malized values of the coupling constants. Then we can 
repeat the integration. In the limit this multi-step pro- 
cedure gives the renormalization group equations for the 
coupling constants, as Kosterlitz suggested. On the dia- 
grammatic language the procedure means that we gradu- 
ally substitute blocks corresponding to small separations 
between the particles by their effective values relative to 
larger scales. The procedure can also be considered as 
increasing an effective size of the defects a — > r. Then 
the renormalization of the coupling constants can be de- 
scribed in terms of the differential renormalization group 
equations. 

At each step of the procedure we deal with correla- 
tion functions like ( |2.9| ). For an interval of scales where 
a variation of the Coulomb consta nt /? is small one can 
use for the function the expression (2.23) where now one 
should substitute the renormalized value of the Coulomb 
constant /3. Analogously, th e re normalized annihilation 
rate is determined by Eq. (3.3) where one should sub- 
stitute the expression (2.21) with the renormalized value 
of the Coulomb constant /3. Next, for the renormalized 
creation rate we should use the relation 



(3.19) 



= / d raS{Ti- T2,r,ra)Rr 



where Ti > T2. The relation ( p.l9| ) can be derived from 
Eq. (|3.6|) if to use the property of S analogous to Eq. 



(3.1). For the renormalized creation rate the relation 
( 3.19 ) is correct only if /3(ri) weakly differs from /3{r2) 
for characteristic values of the parameters ri ~ y/Ti and 
r2 ~ ^M- 



A relation analogous to Eq. ( 3.1£ ) can be formulated 
for the annihilation rate R. The relation leads to the 



same expression (3.3) where /? is now scale-dependent. 
Therefore the renormalized quantity of the annihilation 
constant Ar flows together with /? in accordance with Eq. 
(^.22 ). That is accounted for the non-logari thmic char- 



acter of the integrals leading to the relation (2.22). 

Then we should define the renormalized fugacity in dy- 
namics. For the purpose let us generalize the expression 



dr Rr{r) 



y 



Then Eq. ( 3.16 ) is rewritten as 



y 



(?) 



yo ■ 



(3.20) 



(3.21) 



The relation can be treated as an elementary step of the 
renor ma lization group procedure which is described by 
Eq. (1.5 for y. Thus the renormalization group equation 
for P in dynamics co incid es with one in statics. 

The expression ( 3.17 ) for the correction to the 
Coulomb constant /? can be considered as arising at the 
elementary step of the renormalization group procedure. 
The corresponding renormalization group equation can 
be found if to pass to the differential form and to substi- 
tute ^0 by the renormalized value y in accordance with 
Eq. 



(3.21 



Then we obtain the renormalization group 
equat ion coinciding with Eq. (L5) for /3. The expression 
( p^ for the correction to the diffusion coefhcient leads 



to the following renormalization group equation 



dD 



dln(r/a) 



-y" 



(3.22) 



analogous to the equation (1.5) for /3. We conclude that 
the correction to D is small due to the small value of the 
fugacity and is therefore irrelevant. To avoid a misun- 
derstanding, remind that the variation of the Coulomb 
constant f3 with increasi ng sc ale is also small. Neverthe- 
less, as seen from Eq. (1.5), it is the difference j3 — 2 
that enters the renormalization group equations and the 
variation of the difference can be essential. 



4. CORRELATION FUNCTIONS 

Here, we treat non-simultaneous correlation functions 
of the charge density p (|l.2|) 



F2n{tl, ■ ■ ■ , t2n', fl, . . . , r2„) 
= {p{ti,ri) . . . p{t2n,r2n)) ■ 



(4.1) 
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Note an obvious consequence of the constraint (1.3) 



/ 



(4.2) 



To examine the co rrelation functions ( iA ) we use the rep- 
resentation (1.23). We will assume that all the diagram- 
matic blocks corresponding to small defect-antidefect 
pairs are already included into the renormalization of the 
corresponding coupling constants as discussed in Section 
^. Therefore the fugacity y and the Coulomb constant 
f3 entering all subsequent expressions should be taken at 
the current scale. First we will examine contributions to 
the functions associated with a single defect-antidefect 
pair and then we will consider contributions related to a 
number of defect-antidefect pairs. 



A. Pair Correlation Function 



We start with the pair correlation function 



F2{t2 - <i, r2 - ri) = {p{t2, r2)p{ti,ri)) , 



(4.3) 



with t2 > ti. The average ( |l.3|) can be calculated in ac- 
cordance with the relation ( |1.3C| ) where one should sub- 
stitute the expression ( 1.23| ). We assume U = t2- 

The contribution to the average (4.3) related to a sin- 



gle defect-antidefect pair can be represented by a series 
of the diagrams with two lines constituted of the defect 
and antidefect propagators. The lines start from the cre- 
ation point. A half of the diagrams have the structure 
depicted in Fig. ^ Here we omitted lines and vertices 
corresponding to the interaction of the defects (which is 
implied) and keep only trajectories of the defects. The 
trajectories should pass through the points ri and r2 at 
the time moments ti and t2 (the events are designated 
by black circles). An additional contribution to the aver- 
age (4.3) is determined by similar diagrams where both 
events ti,ri and t2,r2 belong to the same trajectory. 

As above, we dissect the diagram into parts which can 
be treated separately. Let us make a cut along planes 
in r — t space-time corresponding to the time moments 
ti and t2, they are shown in Fig. |9| by dotted li nes. In- 
termediate points appearing in the convolution (3.1) are 
designated in Fig. ^ as and r4. After that the dia- 
gram is divided into two parts separated by the dotted 
line. The part of the diagram to the right from the do tted 
line corresponds to the conditional probability M ( |2.9| ) 
and the part of the diagram to the left from the dotted 
line corresponds to the correlation function 



*2(r-i 



r2 



(V+(t,ri)V'_(t,r2)) 



(4.4) 



The quantity (4.4) can be treated as the probability 
density to find a defect-antidefect pair with a given sep- 
aration. Correspondingly, the integral / d^r$2('") de- 
termines the density of the defect-antidefect pairs. The 



correlation function (4.4) coincides with the integral in 
the left-hand side of Eq. ( 3.20 ). Hence 

<^2{r)^yVr\ (4.5) 

where y is the renormalized fugacity. The equa- 



tions (1.6,4.5) show that asymptotically in the low- 



temperature phase 



^2{r) cx r ^'^ . 



(4.6) 



Note that the same behavior (4.6) is observed up to a 
slowly varying factor in the whole region of scales. 

The diagram depicted in Fig. || gives a convolution 
of <1>2 and M. Adding the contribution corresponding to 
the case where both events ti, ri and ^2, ''2 belong to the 
same trajectory we get the following expression for the 
pair correlation function (4.3) 



F2it,r2-ri) =-2 J d^rg d2r4 $2(^4 - r-i) 

X [M(i,r3,r2,ri,r4)-M(t,r2,r3,ri,r4)] , (4.7) 

where t > 0. At small times t we turn to t he li mit law 
( ^■ll| ). Substituting the expression into Eq. (4/7) we get 



F2{t = 0,ri-r2) 

= 2S{ri - r2) J d^r $2(r-) - 2$2(ri - r2) . (4.8) 



Here the second contribution corresponds to the law (1.7) 
and the term proportional to (5-function is an autocorre- 
lation contribution associated with a single defect. The 
factor in front of the 5-function (which is the density of 
defects) is in accordance with the relation ( |4.2| ). Note 
that at small t the (5-function is converted into a narrow 
bmction of the width ^ \ft. 

It follows from Eqs. ( |2.23| , ^^5p^ ) that for t ~ 



(4.9) 



To justify Eq. (4^) one sh ould check that there are no 
divergences in th e integral (4.7). It can be done directly 
using Eq. ( 2.23 ). The convergence at small separations 
r and vq is accounted for the behavior of the modified 
Bessel functions I,y{x) oc at small values of the argu- 
ment. The convergence at large separations r and Tq can 
be checked using the asymptotic law (2.24). If |f| 3> 
then 



y^{r) 



(4.10) 



r4-2/3|t|/3 • 

The asymptotic law is established in Appendix | 

The behavior of t he pair correlation function deter- 
mined by the laws ( ig| , 4.10 ) corresponds to a conven- 
tional critical dynamics (see, e.g., Ref. |Q) with the dy- 
namical critical index z — 2. However, as we will see 
below, the behavior of the high-order correlation func- 
tions is beyond the conventional scheme. Besides, the 



scaling law i ~ r 
functions as well. 



is true for the high-order correlation 
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B. High-Order Correlation functions 

Here we extend the procedure of the preceding subsec- 
tion to t he case of the high-order correlation functions 
F2n (4J). We wiU assume that h < t2 < ■ ■ ■ < t2n- 
Again, we examine the contribution to associated 
with a single defect-antidefect pair. Corresponding di- 
agrams contain two trajectories starting anywhere and 
passing through the points ri, . . . , r2„ at the time mo- 
ments ti, . . . , t2n- We will designate the trajectories of 
the defect and of the antidefect as x{t) and z{t). 

Let us dissect the diagrams along planes in the r — t 
space-time corresponding to the time moments ti , . . . , 
t2n- Then the diagram is divided into a number of blocks, 
The left block in Fig. 



see Fig. 10 



10 corresponds to the 



object (4.4) and all the other blocks corres pond to the 
correlation function ( |2.9| ). Again, using Eq. {'3A) we can 
write the contribution to the correlation function (4.1) 



associated with a single defect-antidefect pair as the fol- 
lowing convolution 



F2n{h: ■ ■ ■ ,t2n;ri, . . .r2n) 

= n / ^'^J ^'^J *2(a;i - zi) [5{rj 

3 = 1-' 



i)-5{rj-Zj)] 
(4.11) 



Here, one must repl ace th e last factor M{t2n+i ~ ^i-r^i by 



unity. The relation ([4.1l| ) is a generalization of Eq. ([l.7|). 
Thus we got an expression for the correlation function 
which is a multiple integral of functions determined by 
explicit formulas. A recurrent procedure for calculating 
F2n is suggested in Appendi x [c[ 

Of course the expression (4.7 ) for the pair correlation 
function is reproduced by Eq. ( 4.11 ). Note also the fol- 
lowing expression 



lyp(t, ri)p{t, r2)p(0, r3)p(0, Ti)) = 2 [M{t, ri, r2, rs, r^) + M{t, ri, r2, r4, r^)] $2(r3 - T4) 
-B{ri - r2)5{r3 - r^) - B{r3 - r4)6{ri - r2) + J d^r B{r)5{ri - r2)S{r3 - r4) , 



(4.12) 



where 



B(r-3 - r4) = 2 J dVi [M{t, n, r2, rg, r4) 
+M{t,ri,r2,r4,r3)] $2(^-3 - ^4) ■ 



The formula can be found from Eq. (4.11) using the 
relation ( 2.11 ). Naturally, the expression ( 4.12 ) is sym- 
metric under the permutation ri, r2_<-> 1^3, r^, which can 
be checked using the expressions (|T2|,U). The formula 
(4.12) is in agreement with the general property (4.2). 

The recurrent procedure for calculating suggested 
in Appendix ^ shows that there are no divergences in 
the integrals at all steps of calculating F2„. This means 
that we can evaluate the correlation functions from naive 
dimension estimates. Namely, if all space separations 
among jr^ — rj\ are of the same order r* and all time 
intervals are of the order r'i then 



F2n ~ y^{r*)rl 



(4.13) 



In the large-scale limi t w hen (3 is saturated we have in 
accordance with Eq. (1.6) 



F2n oc 



4(n-l)-2/3 



(4.14) 



If some space separations among jr^ — rj\ and/or some 
time intervals differ strongly then one can fo rmulate some 
simple rules following from Eqs. ( ^.23 4.11 ). Let us give 
some examples. If one of the time intervals r is much 
larger than all values of the squared separations jr.; — p 
then the correlation function behaves lik e F2 n oc t^^ . It 
can be proved like it is done in Appendix A 3 for the pair 
correlation function. For small r there appear contribu- 
tions to F2n short-correlated in space (on scales ~ \/t). 



In the limit r — > the contributions turn into (5-functions, 
as it was for the pair correlation function, see Eq. ([1.^), 
representing an autocorrelation of single defects. If the 
points rj can be divided into two "clouds" with a sepa- 
ration r between the clouds much larger than their sizes 
(and all time intervals are much smaller than ) then the 
principal r-dependence of the corre lati on function i^2n is 
the same as in the function <I'2(^) (fl-5|). 

Remember that the charge density p is related to the 
curl of the gradient of the hexatic angle (p for hexatics 
and to the curl of the gradient of the order parameter 
phase for superfluid films. It is instructive to re-express 
our results in terms of the phase gradient circulations of 
Wip over a closed loop C 



r(t, C) = ^ dr Vv? = 27r J 



'p{t,r) 



(4.15) 



where the second integral is taken over the area inside the 
loop. Correlation functions of F's can be rewritten as in- 
tegrals of the correlation functions i^2n ■ As an example, 
consider the following average 



^'2„ = (r(t,c)r(t + r, c). 

xT[t + {2n - 1)t, C]) . 



(4.16) 



Suppose that the characteristic size of the loop r is large 
enough so that we can assume that /3 is saturated, and 
that T ~ r^. Then the following scaling law is satisfied: 
if r ^ Xr and r X'^t then 



2n 



X 



4-2/3 



2n 



(4.17) 



where X is an arbitrary factor. The law (4.17) is a con- 
sequence of Eq. ( 4.14 ) . It has two striking peculiarities. 



14 



First, it possesses a clear critical dependence. Second, it the estimate (1.8) for the simultaneous correlation func - 
is independent of the order n. 

The procedure described above can be generalized to 
include the external pote ntial . Then one should tak e the 
solution of the equation (2.25) for the functio n ( ] 2.9| ) and 
the corresponding expression for the object ( [l.4[ )! 



C. Many-Pair Contributions 

We have established the contributions to the charge 
density correlation functions F2n associated with a sin- 
gle defect-antidefect pair. Now we are going to discuss 
other contributions to the correlation functions related 
to an arbitrary number of defect-antidefect pairs. Cor- 
respondingly, we should take diagrams with a number of 
trajectories passing through the points ri, r2„ at 
the time moments ti, . . . , t2n- The picture illustrating 
the situation is drawn in Fig. |ll| where black circles cor- 
respond to the arguments of F2„: (ii, ri), . . . , {t2n, f"2n)- 
There we omitted blocks related to short-living defect- 
antidefect pairs regarding that the blocks are already in- 
cluded into the renormalization of the Coulomb constant 

As previously, we can dissect the diagrams along the 
planes in the r — t space-time corresponding to the time 
moments ti, t2n- Then any diagram will be di- 
vided into a number of strips, see Fig. |ll|. The part 
of the diagram within each strip can be treated as the 
corresponding matrix element of the evolution operator 
exp(— J dtH) . Then the contribution to i^2n will be writ- 
ten like ( [l.ll| ) as a convolution of the matrix elements. 
Generally, the matrix elements can be estimated like the 
function (^^): each trajectory segment gives a factor 
which scales as 1/t and t scales as r^. But there are ob- 
vious exceptions from the rule. Namely, going back in 
time we will come to a moment where a defect-antidefect 
pair was created. Again, when we consider small sep- 
arations between the defect and the antidefect (which 
is correct for time moments close to the creation time) 
we can take into account only the Coulomb interaction 
between the two created defects. The corresponding re- 
gions in Fig. |ll| are inside the ovals. Each such region 
produces the factor y^. Therefore generally F2„ oc y'^'', 
where k is the number of the pairs. Taking into account 
also a scale-dependent factor we get 



F2n - y^'ir.y 



(4.18) 



where we assume that all space separations are of the 
order and all time intervals are of the order 



The expression ( |4.1§| ) is a generalization of Eq. ( 4.13 ). 
Comparing these tw o exp ressions we see that t he ra tio 
of the contribution ( [1.18| ) to the contribution ( 4.15 ) is 
the (fc — l)-th power of a dimensionless small parameter 
2/^(r*). Thus we conclude that the leading contribution 
to F2n is related to a single defect-antidefect pair, that 
corresponds to fc = 1. Now we can explain the origin of 



tions which obviously does not coincide with Eq. ( 4.13 ). 
The estimate (1.8) in terms of Eq. (4.18) corresponds to 
k = n. The reason is quite obvious: Two defects can- 
not pass simultaneously through 2n points and at least 
k = n defect-antidefect pairs should be taken to obtain 
a non-zero contribution to the simultaneous correlation 
function i^2n- The situation is illustrated by Fig. ||. Note 
that the estimate ( p,.^ ) is not correct for the autocorrela- 
tion contributions proportional to ^-functions, as written 
in Eq. (Q. 

Thus we have two different regimes: for simultaneous 
and for non-simultaneous correlation functions. Let us 
establish the boundary between the regimes. For the pur- 
pose we should consider small time intervals where the 
single-pair contribution is finite but small. The small- 
ness is associated wi th di ffusive exponents presented, e.g., 
in the expression (2.23). Therefore the characteristic 
time where the simultaneous regime passes into the non- 
simultaneous one can be estimated as 



t ' 



I ln[2/(r) 



(4.19) 



where r is a space separation corresponding to the small 
time interval. In the low-temperature phase on large 
scales (where (5 is saturated) we have |lny| sa (/? — 
2) ln(r/a). 



5. SUPERFLUID FILMS 

Let us consider superfluid films. The equation of mo- 
tion for the vortices contains an ad ditional term (Magnus 
force). Thus instead of Eq. ( 1.10 ) we should write (see 
Ref. 0) 



D 

T 



(5.1) 



where 7 is a new dimensionless parameter. The equation 
( ^.iD can be derived in the spirit of the procedure pro- 
posed by Hall and Vinen for the id superfluid, see Ref. 
p5[ . Huber |2^ argued that the same equation is correct 
for spin vortices in planar 2d magnetics. 



For the superfluid films the "charge density" (^_^) is 
proportional to the vo rtici ty curlDg. To calculate cor- 
relation functions F2„ (4T) one can use the scheme de- 
veloped in the previous section s. T he only difference is 



that instead of the expression ( 2.23 ) for the conditional 
probability M (2.9) one should use the solution of the 
equation 



-Rir)M + S{t)S{r -rQ)S (q 



(5.2) 



The variables r, ro and g are introduced by Eq. ( ^.12 ). 
Again, on scales r ^ a one can omit the term with the 
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annihilation rate R in Eq. (5.2) demanding a finite value 
of M at r — ^ instead. 

Unfortunately, a cross term over r and g app ears in 
the operator in the right-hand side of Eq. (5.2). Thus 
one cannot o btain an explicit expression for M of the 
type of Eq. (2.23). Nevertheless, this additional term 
has the same dimensionality as the other terms and does 



not change the scaling estimates M ~ < , t ' 



deter- 



mining the function M . Moreover, the equation for the 
object 



Sit,r,rn)= J g M{t,r, Q,ro,r3/2 + r4/2) 



(5.3) 



following from (|5.2| ) is identical to Eq. (2.17). The refore 
for the object ( ^.3D w e have the same series ( 2.16 ) with 



the coefficients (2.21). 

Looking through the derivation pre sented in Section ^ 
we see that just the function S ( |5.3| ) enters all the re- 
lations. Therefore we can make the same assertions as 
previously. First, on large scales the annihilation coef- 
ficient A is equal to its universal v alue (2.22). Second, 



we c an write the same expression (4.5) for the average 
(4.4). Third, in dynamics we get the same renormaliza- 
tion group equation (1.5) for /3, see Appendix ^. Fourth, 
the renormalization of the diffusion coefRcient is irrele- 
vant. And finally, one can assert that a re normalization 
of the parameter 7 introduced by Eq. (5.1) is determined 
by the equation 

d7 2 



d ln(r/a) 



analogous to ( p. 22 ). Therefore the renormalization of 7 
is irrelevant. Again, the scheme can be generalized to in- 
clude the "external potential" , which now is the average 
value of the superfluid velocity. 



Next, we proceed to the correlation functions F2„ (4.1 ). 
Formally they are determined by the same convolution 
(4.11) as previously. However, on e s hould substitute 
there the solution of the equation (^.2|). Therefore the 
concrete expressions for will be di fferent. Neverthe- 
less, the estimates like ( 4.13 , 4.14|4.18| ) remains true be- 
cause of the following reasons. First, due to the same 
dimensionality of all the terms in the right-hand side of 
Eq. ( ^ ) the function M possesses the simple scahng 
properties noted above. S econ d, there are no divergences 
in the convolutions like ( 4.11 ) determining the objects. 
To prove the second property, wc should analyze a be- 
havior of M at large and at small separations. In the 
case rrp /t ^ 1 the characteristic values of the separations 
ri— rs and of r2— r4 are ^ \ft and are consequently much 
smaller than \ri — (or — r4|). Then it is possible 
to neglect all terms containing |ri — r2| in denominators 
in the equation ( ^.2[ ) and we come t o a p urely diffusive 
equation leading to the asymptotics (2.24). It is possible 
to establish that the small-scale of the conditional prob- 
ability M for the vortices coincides with that examined 



above. The properties ensure convergence of all interme- 
diate integrals appearing at calculating the correlation 
functions of vorticity 



We have also t he sam e scaling law ( 4.17 ) for the cor- 
relation function (4.16) of the integrals (4.15) which are 
now proportional to the circulations of the superfluid ve- 
locity. We conclude that all the scaling laws for the cor- 
relation functions of the vorticity and their asymptotic 
behavior remains the same as previously. 



6. DISCUSSION 

The main result of our consideration is the expression 
(4.13) for high -order correlation functions of the "charge 
density" (4.1) which is disclinicity for hexatic film s and 
vorticity for superfluid films. We see from Eq. (4.13) that 
the high-order correlation functions arc much larger than 
their normal estimates via the pai r correla tion function. 
Namely, in accordance with Eqs. (O 4.13| ) we have 

F2„/F2" ^ » 1 , (6.1) 

where y is the renormalized fugacity. The asymptotic 
behavior of the ratio at large scales is determined by the 
law ( |l.6| ). Though at developing our scheme we accepted 
that the defect-antidefect pairs constitute a dilute solu- 
tion we hope that the scaling law (6.1) is universal. The 
ground for the hope is the renormalization group proce- 
dure (formulated in Ref. ||l^) which shows that on large 
scales we come to an effectively dilute solution of the 
pairs. We believe that the most interesting fact to be 
compared with experiment or numerics is the scaling law 
(4.17) which is a consequence of Eq. (|6.1|) . 

The physics behind the inequality ( |6.1| ) is as follows. 
The main contribution to the correlation functions is as- 
sociated with a single defect-antidefect pair. Though 
the contribution associated with a number of defect- 
antidefect pairs contains an additional huge entropy fac- 
tor it has also an additional small factor associated with 
small probability to observe defect-antidefect pairs with 
separations larger than the core radius a. The smallness 
is accounted for the strong Coulomb attraction. The con- 
sidered effect is a consequence of the competition of those 
two factors. The result of the competition manifests in 



the law (4.18) which gives the estimate for the contri- 
bution associated with k defect-antidefect pairs. For si- 
multaneous correlation functions nothi ng s imilar occurs 
and we have the conventional estimate (1.8). That is the 
reason why the effect cannot be observed in statics. The 
property is directly related to causality since a defect- 
antidefect pair cannot simultaneously pass through 2n 
points and at least n defect-antidefect pairs is needed to 
get a non-zero contribution to the simultaneous correla- 
tion function F2„, see Fig. That explains the esti- 
mate (1.8). Thus we have two different regimes for si- 
multaneous and non-simultaneous correlation functions. 
The characteristic bound ary ti me separating those two 
regimes is written in Eq. (4.19). 
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The considered effect resembles intermittency in tur- 
bulence (see, e.g., Ref. Il4[) le ading to large r-dependent 



factors in the ratios like (6.1) in the velocity correlation 
functions of a turbulent flow. However, as is seen from 
Eq. (6.1), for the defects the large r-dependent factors 



are related to the ultraviolet cutoff parameter a whereas 
for intermittency in turbulence the large r-dependent fac- 
tors are related to the infrared (pumping) scale. Our sit- 



uation is thus closer to the inverse cascade (see Ref. |27 



realized on scales much larger than the pumping length. 
There are experimental data [ p8| concerning the inverse 
cascade in 2d hydrodynamics and analytical observations 
concerning the inverse cascade for a compressible fluid 
pof which indicate the absence of the intermittency in 
the inverse cascades. Note that only simultaneous ob- 
jects were examined in the works p8| , p9t , and there is 
no intermittency in our simultaneous correlation func- 
tions. So, based on the analogy, one may think that 
for the inverse cascades non-simultaneous objects reveal 
some intermittency. 

The consideration presented in our work is applicable 
to superfluid films. There exist also films and quasi-2c? 
systems of different symmetry. Huber p^ ] argued that 
the same equation as for quantum vortices is correct for 
spin vortices in planar 2d magnet s. W e believe that our 
approach based on the equation ( I.IOD is correct for the 
dynamics of disclinations in hexatic films like membranes, 
freely suspended films and Langmuir films. Next, the 
above scheme seems to work also for dislocations in solid 
films. The system needs a special treatment since a mod- 
ification should be introduced into the procedure. Maybe 
some features of the presented picture can be observed 
also in superconductive materials, especially in high-Tc 
superconductors. There are analytical and numerical in- 
dications that for a purely Langevin dynamics of the or- 
der param eter there are logarithmic corrections to the 
law ( pTTo| ), see Refs. |o|-|§. We believe that the loga- 
rithms are destroyed if to switch on an interaction of the 
order parameter with other degrees of freedom. Never- 
theless, it would be interesting to generalize our scheme 
including the logarithmic corrections. 
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APPENDIX A: 

In this Appendix we present some calculations leading 
to the results presented in the main body. 



1. Some relations for the pair conditional probability 

Here, we dedu ce s ome relations concerning the corre- 
lation function (2^). We examine the solution of the 
equations (2.17) with the annihilation term R omitted. 
Then we should accept a finite value of Sm at r ^ as 
is explained in the main text. 

Performing the Fourier transform over t we get from 
Eq. (pT7|) 



-'S'm(w) 



1 Til 

a2 + (l + 2/3)-<9,--^ 



+ 



d{r - ro) . 



(Al) 



The solution of the equation (Al) which tends to zero at 
r — > cxD and remains finite at r — > can be written as 




if r > ro , 



(A2) 



where = ^JfW^^^n? and we used the relation 

K,{z)dJ,{z) - I,{z)d,K,{z) = . 

At the next step one performs the inverse Fourier trans- 
form. Deforming the integration contour to the negative 
imaginary semi-axis and using the relation 6.633.2 from 
Ref. 13 one gets Eq. (|^). 

Next, we present two integral relations for the condi- 
tional probability (2.9) in the above approximation. Let 
us calculate the total probability to find a defect and 
an antidefect in any points at a fixed time separation t. 
Using the relation 6.631.1 from Ref. [B^ we get 



j dVi dVaAf = j A^rS{t,r,ro) 



(A3) 



r(i + /3) V8t 



exp 



8t 



1^1 + 



We conclude from Eq. (A3) that the total probability 
diminishes with increasing time t. It is accounted for 
annihilation of defects at collisions. Usi ng th e relation 
6.611.4 written in Ref. we find from (p!23|) 
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I 



X ■ 



dtM 



1 



m— — OO 



+ 00 



{2rror 



^^2 „ 4,.2^2 1^^- ^ ^^2 _ 4,-2,. 2^ 



(A4) 



Tq and V 



where = (j"! + ^2 — '"a — + 

+ m^. The expression determines a distribution of 
the defect and of the antidefect over space provided the 
pair was created near the origin at any time. 

No w we calculate the integral (3.5). Substituting there 
Eqs. (3.3, A3|) we get after integrating in part 



-stt y dz 




r(i + /3) 



iFi(l,/3+l;z)-l 



(A5) 



where z = rg/Sr. One can easily check using the asymp- 
totic exp ress ion for iFi(l,/? + l;z) at large z that the 
integ ral (A5) converges. Let us now rewrite the expres- 
sion (A5) as 



Ar = — Stt lim / dz 



r(i + /3) 



exp(— z — az) 



Xii^i(l,/3 -I- 1; z) — exp(— az) 



(A6) 



Now integrals of both contributions to the integrand can 
be found explicitly (see the relation 7.621.5 from Ref. 
p7t). P assin g then to the limit a ^ we come to the 
answer ( 2.22 ). 



2. Asymptotics of the Pair Conditional Probability 

Here, we will be interested in the asymptotic behavior 
of the expression ( 2.23| ) where the parameter z — rro/(4i) 
is large; z 3> 1. More precisely, we will examine the func- 
tion 



-t-00 



(^) = ^ exp(imv5) (z) , 



(A7) 



It is obvious that T{z,ip) — 



where v 

T{z, —(p). Below we will assume ip > 0. 

First of all, we convert the sum over m in Eq. (A?) 
into an integral using the identity 

+00 ^ „ 

f{m) — 7: dm coth(— i7rTO)/(m) 
^ / 1 

m — — 00 ~r 



— — J dm coth(— i7rm)/(m) , 

where the contour in the first integral goes above the real 
axis and the contour in the second integral goes below the 



real axis. One can easily check that if /(— m) is complex 
conjugated to /(m) then the second integral is equal to 
the complex conjugated of the first integral with the sign 
minus. Therefore we can write 

r(z. If) = Re J dm coth(— iTrm) exp{imtp) ly (z) . (A8) 

Next, we are going to shift the integration contour in 
Eq. (A8) into the upper m-semi-plane. Then we should 
separately consider two contributions to ly which are de- 
termined by two terms in its integral representation 



1 r 

Iv{x) = — dS exp(x 

Jo 



cosi9) cos(i^i9) , 



(A9) 



Substit utin g the first contribution in the right-hand side 
of Eq. (|A9|) we get 



Ti{z,(p) = Re J dm coth(— iTrm) 

d-d 

— exp(z cosi}) exp{imip) cos^vd) . 



(AlO) 



As we will see, the answer will be determined by a vicinity 
of a saddle point where —im 3> 1. There one can substi- 
tute coth(— i7rm) by unity and v by m. Then one obtains 
from Eq. ( |A10| ) the following saddle-point conditions 



i9 



iz sin i9 . 



One can easily find from Eq. ( AlO ) in the saddle-point 
approximation 



Ti {z,ip) « exp(z cos Lp) . 



(All) 



where we have taken into account an pre-exponent be- 
sides the exponent. The above scheme is obviously bro- 
ken at small </? or at <p close to tt since the saddle-point 
value of m tends to zero there. Thus the cases need a 
separate consideration. 

Let us consider the case wher e ip is close to zero. Then 
the main contribution to Eq. ( [A10| ) is gained from the 
region of integration over near zero. There it is possi- 
ble to expand the factor at z as cosd « 1 — 19^/2. Then 
the integration over Lp is performed explicitly. Returning 
also to the sum we get 



Ti(z,(/?) 



exp z 



+00 

E 



exp imip 



m 
"27 



where we substituted exp(— /3^/2z) by unity. Since both 
(y5 <C 1 and 1/z <C 1 we can replace here summation by 
the integration over m. After the substitution we get 

ri(z, ip) w exp (z — zkp^ jlL) , 



what reproduces Eq. (All) at small i-p. Analogously the 
case of (y9 close to tt can be examined. Again, we repro- 
duce Eq. (|aTT|). 
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Next, we examine the contribution to the function T 



associated with the second term in (A9). For the purpose 
we return to the initial expression ([A7|). The second con- 
tribution can be written as 

T2(z, V') — / dt exp(im(y9 — vt) 

^ m— — 00 

X sm{'Kv) exp(— zcosht) . 

If the sum is dominated by first terms in the sum over m 
then we have 



\T2{z,ip)\ < ^exp(-z) , 



(A12) 



3. Asymptotics of the pair correlation function 

Let us consider the asymptotics of F2 for times t much 
larger than rf2- Then th e ch aracteristic values of both 
r and ro in the integral (4/7) are much larger than ri2 
and therefore we can neglect a dependen ce on ri2 t here . 
Substituting the expressions (UUItI) into Eq. Q 
and then neglecting ri2 we get after integrating over the 
angle between r and rp 



F2 



dr dro 

y,4-2/3;2 / J./3-1 0-1 



y 



exp 



At 



00 



which can be obtained after substituting cosht w 
into the integral. Thus we should estimate the contribu- 
tion determined by large |to|. Then 



sin(7ri^) sa (— l)*^ 



2|m| 



Substituting the expression into the sum we get a contri- 
bution to T2(z, ip) 



I dt exp(— zcosht) 



0. 

E- — - cos[(7r — ip)m] exp(—mt) 



/•OG 

= / di exp(-2Cosht) $ [t^ -f (tt - (^)^] , 
Jo 

where we introduced a function which behaves as 

<i>(x) ~ Inx if a; < 1 and |<i>| < 1 if x > 1. Therefore 
the contribution determined by large |m| reproduces the 
same estimate (A12). The est imat e ( A12| ) shows that 
\T2\ is much sm aller than Eq. (All). We conclude that 
the expressi on ( All ) determines the main contribution 
to the sum (A7). 



Rewriting the answer (All) for the sum (A7) in terms 
of the original variables we get 



-t-CX3 



exp(im(^) I„ w exp (J^) ■ (A13) 



m— — cxD 



Then the expression (A13) leads to 



M 



exp 



ro)' 



8t 



X exp ■ 



(ri + r2 - r3 - Vj) 
8t 



(A14) 



where 1/ — (3^ + (2m -|- 1)^ and tq = ri — r^, r = 
— r2- The ratio /r'^~^l^ is a slowly varying function 
and therefore it is placed outside the integral. Note that 
the integral is nonzero due to Ii, > 0. Performing here 
the integration over r/rp we get 



F2 



la V 

00 

l2m+i{y)lu{y) 



m— — 00 



where we used the relation 3.547.4 from Ref. ||37||. The 
integral over y converges at any m obviously. Thus we 
should examine the convergence of the sum over m. At 
large v we can use the asymptotics 



X exp 



1 



1 



/2^ (a;2 + j.2)i/4 



(A15) 



and the integral over y is gained at large y, where the 
asymptotics Ko{2y) w y^n/Ay exp(— 2j/) works. There- 
fore 

^"-^ifo(2y)/2„.+i {y)iAy) 

■■ exp(2i//) , 



Jo z/3-i/2Vi^ 



f{z) = - z + In 



1 + Vz^ + 1 



Since the function f{z) increases monotonically, the in- 
tegral over z is gained where f ^ that is at z ~ v. 
That leads to an estimate 



V 



dy 

—Ko{2y)l2m+i (y) lu (y) 



1 



7,2/3-1 



The estimate means that the sum over m converges at 
l arge m since /3 > 2. Thus we come to the estimate 

(H). 
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APPENDIX B: CORRECTION TO THE 
COULOMB COUPLING CONSTANT 

Here, we calculate the correction to the Coulomb cou- 
pling constant P produced by small defect-antidefect 
pairs. We start from Eqs. ( 3.1C , 3.11 ). The correction 
can be found from the relation 



exp 



+ 



dtAH2(0,CT) 



(Bl) 



exp 



dtd^r 



where the angular brackets mean averaging over the small 
defect-antidefect pairs and (f), a in the exponent are to 
be treated as "external" fields. Really, one of the con- 
tributions to ATC is depicted in Fig. ^ where the upper 
arrowed dashed line correspond to the field (j) and the 
lower arrowed dashed line correspond to the field a. The 
other contribution to ATi. is associated with a similar di- 
agram where both "external" dashed line are attached 
to the same trajectory. It is explained in Sec. |^ what 
objects correspond to the different parts of the diagram. 
Here we present an analytical expression written in ac- 
cordance with the two noted diagrams (see designations 
of the points in Fig. ||) 

= 2 y dt2 dta (Pxi (fx2 d^xg (fx^ 

X Vla*2(a;i - X2)M{t3 - t2,X3, X4, Xi,X2) , (B2) 



where the objects $2 and M are defined by Eqs. (2^,4.4). 
The explicit expressions for the functions are written in 
Eqs. ( ^!23| , p6| ). 

We are interested in the situation when the "external" 
fields (f) and a in Eq. (B2) vary on scales much larger 



than a characteristic size of the defect-antidefect pair. 
Then one can substitute 

a{t3,X4) - (7{tz,Xz) V^(j{t2,Xi)[Xi^ - Xz-^) , 

where we expanded the difference and then shift the ar- 
gument of the factor Vcr. Substituting the expression 
into Eq. (B2) and integrating over q = {x^ 4- 0:4 )/2 we 
get 

AH2 = ^'^ j dtd^rVa(t){t,r)\7^a{t,r) (B3) 
X J dTd^xd^roXjVa'^2{ro)S{T,x,ro) , 

where Vq = Xi — X2, x = X3 — X4, t = t2, r — ^3 — t2, 
r — Xi, and we substituted the expression (2.13). Next, 
the integral ( [B3| ) is proportional to Sap (due to averag- 
i ng over an gles). Taking into account also the expressions 
( pjjpl^ ) we get 



An2=iTT^f3 J dtd^r\/(j){t,r)S/a{t,r) 

00 

X y dTdxdrox'^<^2{ro)S{T,x,ro) . (B4) 



Comparing the expression with Eq. ( ^.11| ) we conclude 
that the integral in Eq. (B4) can be considered as a 
correction to the coefficient l/(47r/3). Thus we get 



Af3 = -IGir^P^ J dTdxdrox^^2iro)S{T,x,ro) . (B5) 




Substituting here the expression ( ^.21 ) we obtain 



A/3: 



27r2/33 



dx dro a;^$2(ro) 



{x^ + rl - \x^-rl\y (r^^H 
{2xroY 



(B6) 



where v — + \ and we used the relation 6 .623. 3 

from Ref . |Q . Now, remembering the expression ( j3.16D , 
one can easily check that the integral over rp in Eq. (pq) 
converges both for small and large tq. Therefore t he in - 
tegral is gained at tq ~ a;. Thus, substituting Eq. (3.16), 
we get from Eq. (B6) the answer ( |3.17| ). 

Now some words about generalization to the case of 
superfluid films discussed in Sec tion ^. The new term in 
the right-hand side of Eq. ( |5.l|) produces an additional 
contributi on to the "Hamiltonian" of the system. Then 
the term (3.10) should be substituted by 



Hi 



dV 



7/;+Va7/>+ - 0_VaV_ j Va0 (B7) 

a 



The expression (B7) means a modification of the triple 
vertices shown by the arrows in diagrams, see e.g. Fig. 
^. Next, we can repeat the same steps as above starting 
from the Hamiltonian determined by Eq. (B7). Then the 
intermediate expressions like Eqs. (B^,B3) will be modi- 
fied. The terms which does not contain 7 are reduced to 
the c ond itional probability S and give the same correc- 
tion (B6) to /?. The terms containing 7 do not produce 
relevant terms at all because of the identity 



d^r VaCrea/3V/30 = . 



Thus for the superfluid films we get the same answers 

HtI). 
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APPENDIX C: RECURRENT PROCEDURE 



We see from the expression (4.11) that the correlation 
function i^2ra can be calculated step by step: first inte- 
grating over Xi,zi, then over X2,Z2 and so further. To 
formulate the procedure one should introduce auxiliary 
objects 

Hk(ti,...,tk;ri,.. . ,rfc_i;a;fe,2fc) 

fe-i 



d^Xj d^Zj [S{rj — Xj) — S{rj 



xM{tk - tk^i,Xk,Zk,Xk^i,Zk-^i) . . . 

xM{t2 - ti,X2,Z2,Xi,Zi)<p2{Xi - Zi) 



(CI) 

One can easily check that the functions are sym- 
metric under permuting x^ and Zk for odd k and are 
antisymmetric under the permutation for even k. That 
depends on the number of factors containing differences 
of 5-functions which are antisymmetric under permuting 
Xj and Zj whereas the factors M are symmetric under 
the permutation. The correlation fmiction i^2n can be 
restored from H2n'- 



-Pzn = 2 j (l^QH2n{ti]ri,...,r2n-l,r2n;Q), 



(C2) 



as follows from Eq. (4.11). Note that an attempt to cal- 
culate odd correlation functions F2n+i from H2n+i using 
the same scheme gives zero (as it should be) since i?2n+i 
is symmetric under permuting X2n+i and Z2n+i whereas 
the difference S{r-'X2n+i)—S{r—Z2n+i) is antisymmetric 
under the permutation. 

To obtain Hk we can us e a recurrent scheme. First of 
all we conclude from Eq. (CI) that 

Hi{t,x,z)^^2ix-z). (C3) 

Next, one can easily obtain from the definition ( |Cl| ) the 
following recurrent relation 

Hk+i{ti, ■ ■ ■ , ife+i; ri, . . . , rk,Xk+i, Zk+i) 

(fxk (fzk [S{rk - Xk) - 5{rk - Zk)] 

xM{tk+i — tk,Xk+i, Zk+i,Xk, Zk) 

xHkih, . . . ,tk;ri, . . . ,rk^i,Xk,Zk) ■ (C4) 

Taking into account the sym metry properties of we 
can rewrite the relation (C4) as 



H2n{il, ■ ■ ■ , t2n', ^1, • ■ • , ^2,1-1 , X2n , Z2n) — — J d^g i?2n-l (^1 , • ■ • , ^2n-i; '"l, • ■ • , f"2n-l , ■''2n , Q) 

X [M{t2n - t2n~l,X2n,Z2n, Q,r2n-l) - M{t2n ~ t2n~l,X2n,Z2n,r2n-l, Q)] , (3.5) 

i72„+i(ti,...,W;n,...,r2„,a.2„+i,.2„+i) = /d^^.i/2„(ti,...,t2„;n,...,r^^^ 

X [M{t2n+1 - hn,X2n+l, Z2n+1, Q, r2n) + M{t2n+1 - hn, X2n+1, Z2n+l,r2n, Q)] ■ (3.6) 



Thus we can subsequently obtain the functions Hk and 
then get the correlation functions i^2n in accordance with 
Eq. (g). 

Now we are going to establish scaling properties of the 
correlation functions. As is seen from the recurrent rela- 
tions ( 3.5.3^ and Eq. (C2) the crucial point is conver- 
gence of integrals over q there. There are no prob lems 



with finite p since one can conclude from Eq. ( 2.23 ) that 
M remains finite at any val ues o f its space arguments. 
Indeed, it follows from Eq. ( 2.23 ) that the only danger- 
ous case is r ^ 0. Then the arguments of also tend to 
zero. Substituting there the first terms of the expansion 
I^{x) ~ (a;/2)''/r(l + v) we conclude that the value of 
M remains finite at r ^ 0. Therefore w e should examine 
a behavior of the integrands in Eqs. (3.5,3.6) at larg e 
Q. The behavior is dominated by the asymptotics ( 2.24 ). 
First of all note that at large g 



Hk{ti, . . .,tk;ri, . . .,rk-i,Xk, q) cx g 



(3.7) 



The behavior (3.7) is characteristic of $2, that is char- 
acteristic of Hi, see Eqs. (03,3.16). Then the behavior 
i s repro duced at each step of the recurrent pr ocedure 
(3.5 3^). Indeed, in the first terms in Eqs. (3.5,3.(:) 



M ^ 1 if ^ is close to x and M exponentially decays 
at incre asing g — x. Therefore at large x the integrals 
over Q are gained near x, that leads to the above 
assertion. The asympto tic behavior (3.7) ensures con- 
vergence of the integral (C2) determining the correlation 
functions. 

Thus there are no divergences in the int egrals a t all 
step s of calculating i<2n in accordance with (3. 5 3. 6) and 
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FIG. 1. Trajectories passing through given points. 



22 



FIG. 2. Possible and impossible trajectories passing through four points at a given time moment. 




FIG. 3. Typical diagram block. 
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FIG. 4. Typical diagram for the conditional probability M. 
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FIG. 5. Vicinities of Creation and Annihilation Points 



FIG. 6. Typical diagram contributing to renormalization of the effective "dielectric constant" 



FIG. 7. A more complicated diagram giving a correction to the Coulomb interaction. 
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FIG. 8. Illustration to the renormalization of the diffusion coefficient. 
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FIG. 9. Trajectories contributing to the pair correlation function. 
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FIG. 10. Trajectories passing through a number of points. 
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FIG. 11. A number of defect-antidefect pairs passing through given points. 
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